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ABSTRACT 

Context. This paper is part of a larger project in which we systematically study the chemical abundances of extra-galactic post- 
asymptotic giant branch (post-AGB) stars. The aim of our programme is to derive chemical abundances of stars covering a large range 
in luminosity and metallicity with the ultimate goal of testing, constraining and improving our knowledge of the poorly understood 
AGB phase, especially the third dredge-up mixing processes and associated s-process nucleosynthesis. 

Aims. Post-AGB photospheres are dominated by atomic lines and indicate the effects of internal chemical enrichment processes 
over the entire stellar lifetime. In this paper, we study two carefully selected post-AGB stars: J051213.81-693537.1 and J051848.86- 
700246.9 in the Large Magellanic Cloud (LMC). Both objects show signs of s-process enhancement. The combination of favourable 
atmospheric parameters for detailed abundance studies and their known distances (and hence luminosities and initial masses) make 
these objects ideal probes of the AGB third dredge-up and s-process nucleosynthesis in that they provide observational constraints for 
theoretical AGB models. 

Methods. We use high-resolution optical UVES spectra to determine accurate stellar parameters and subsequently perform detailed 
elemental abundance studies of post-AGB stars. Additionally, we use available photometric data covering optical and IR bands to 
construct spectral energy distributions for reddening and luminosity determinations. We then estimate initial masses from theoretical 
post-AGB tracks. 

Results. We obtained accurate atmospheric parameters for J051213.81-693537.1 (T e ff = 5875 + 125 K, log g = 1.00 + 0.25 dex, 
[Fe/H]= -0.56 + 0.16 dex) and J051848.86-700246.9 (T eff = 6000 ± 125 K, log g = 0.50 ± 0.25 dex, [Fe/H]= -1.06 ± 0.17 dex). 
Both stars show extreme s-process enrichment associated with relatively low C/O ratios of 1.26 + 0.40 and 1.29 + 0.30 for J051213- 
693537.1 and J051848-700246.9, respectively. We could only derive upper limits of the lead (Pb) abundance. These upper limits show 
a possible very slight Pb overabundance with respect to heavy 5-elements for J051213-693537.1, while J051848-700246.9 shows an 
upper limit of the Pb abundance similar to [hs/Fe]. A comparison with theoretical post-AGB evolutionary tracks in the HR-diagram 
reveals that both stars have low initial masses between 1.0 and 1.5 M 0 . 

Conclusions. This study adds to the results obtained so far on a very limited number of 5-process enriched stars in the Magellanic 
Clouds. With the addition of the two stars in this study, we find an increasing discrepancy between observed and predicted Pb 
abundances towards lower metallicities of the studied 5-process rich post-AGB stars in the Magellanic Clouds. The more metal-rich 
J051213-693537.1 fits the theoretical Pb abundance predictions well, while the five other objects with [Fe/H] < 1, including J051848- 
700246.9, have much lower Pb overabundances than predicted. In all objects found so far, including the objects in this study, the C/O 
ratio is very moderate because of the enhancement of O as well as C. We find that all 5-process rich stars in the LMC and SMC studied 
so far, cluster in the same region of the HR-diagram and are associated with low-mass stars with a low metallicity on average. We 
corroborate the published lack of correlation between the metallicity and the neutron irradiation, while the neutron exposure ([hs/ls]) 
is strongly correlated with the third dredge-up efficiency ([s/Fe]). These correlations seem to hold in our Galaxy as well as in the 
Magellanic Clouds. 
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1. Introduction 

Post-asymptotic giant branch (post-AGB) stars are low- to 
intermediate-mass stars (M < 8M 0 ) in a transient evolutionary 
phase between the asymptotic giant branch (AGB) and the white 
dwarf (WD) phase. In single stars, a superwind mass loss termi¬ 
nates the AGB phase after which the star evolves onto the post- 
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* Based on observations collected with the Very Large Telescope at 
the ESO Paranal Observatory (Chili) of programme number 088.D- 
0433. 


AGB phase, eventually cooling down as a white dwarf. During 
the AGB phase, the end products of internal chemical processes 
like C,N,0 and 5-process elements are transported to the stellar 
surface by multiple dredge-up events, called the third dredge- 
up (TDU). The 5-process synthesis in AGB stars is an impor¬ 
tant contributor to the cosmic abundances past the iron peak and 
these stars are also expected to be important contributors to the 
cosmic carbon and nitrogen enrichment (e.g Romano e t al.|2010[ 
[Kobayashi et al.|20TT| . 

Photospheric spectra of AGB stars are dominated by molec¬ 
ular bands and they are often veiled by circumstellar dust, mak¬ 
ing it difficult to study abundances of individual elements (e.g. 
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Abia et al.|2008| ). However, post-AGB spectra are dominated by 
atomic transitions, allowing extensive abundance studies of in¬ 
dividual elements. Their photospheres bear witness to the total 
chemical changes accumulated during the stellar lifetime, mak¬ 
ing them ideal tracers for the study of AGB nucleosynthesis and 
associated TDUs. 

Chemical studies of Galactic post-AGB stars s how a large 
variety of chemical abundance patterns (e.g. |Van Winckel|2003[ 
and references therein). Only about 25% of the studied Galactic 
post -AGB stars indicate strong ^-process enrichment (e.g. |Rao| 
et al.|2012j and references therein). Although stars in the Galac¬ 


tic sample of post-AGB stars are well studied, their unknown dis 
tances hamper accurate luminosity determinations and hence ac 
curate initial mass estimates, which are a key parameter in con 
straining single star evolution and nucleosynthesis AGB models. 

One way to overcome this distance problem, is to study post- 
AGB stars in the Magellanic Clouds. The distances to the Mag¬ 
ellanic Clouds are well known, with the LMC at a distance of 
about 50 kpc ( [Keller & Wood|2006{ [Reid & Parker|2010[|Storm 


et al.|2011~]) and the SMC at a distance of about 60 kpc (Keller & 


Wood|2006| . Furthermore, the Magellanic Cloud stars generally 


have lower mean metallicities than their Galactic counterparts, 
with a mean metallicity of ~ -0.5 dex for the LMC (e.g. |Geisler| 
et al.||2009t |Lapenna et al.||2012|) and ~ -0.7 dex for the SMC 


( Luck et al.|1998| ). This allows for the study of post-AGB stars 
in a broad metallicity range. 

Therefore, we initiated a larger project in which we search 
for post-AGB stars in the Magellanic Clouds and then use these 
objects to study the poorly understood AGB TDUs and s-process 
nucleosynthesis. [van Aarle et~ak| ( [2011| ) and |Kamath et al.| ( [2015| 
accepted), provide catalogues of spectroscopically verified op¬ 
tically visible post-AGB candidates in the Large Magellanic 
Cloud (LMC) and Ka math et al.| ( |2014| ) provides a catalogue of 
spectroscopically verified optically visible post-AGB candidates 
in the Small Magellanic Cloud (SMC). 

In this paper, we present the detailed chemical study of two 
newly identified LMC post-AGB candidates that were carefully 
selected (based on their stellar parameters and spectral features) 
from our optically visible post-AGB LMC catalogues. 

Our previous abundance studies of strongly ^-process en¬ 
riched post-AGB stars in the Magellani c Clouds by |De Sm edt 


[et al.| ( [Ml2| ), lvan Aarle et al.| ( |20131 ), and [De Smedt et al.| ( |2014| ) 
show that there are still large discrepancies between observed 
abundances and theoretical predictions. The biggest discrepan¬ 
cies are related to the C/O ratio since observationally derived 
C/O ratios are small in comparison to the theoretical predic¬ 
tions. The theoretical carbon abundance is often overestimated, 
whereas the derived oxygen abundance is larger than predicted. 
Also lead (Pb, Z=82), which is considered to be the end product 
of the ^-process chain due to its double magic number, is strongly 


overestimated by theoretical models as shown in De Smedt et al.| 
|20l4l ). 

We describe the observations and data reduction of both ob¬ 
jects in Sect. [2] Detailed abundance analysis of the obtained 
spectra is presented in Sect. [3] followed by the discussion of 
the abundances in Sect. [4] Luminosity determinations and ini¬ 
tial mass estimates are presented in [5 We discuss the neutron 
irradiation of our sample stars in Sect. 6] Finally, we end with a 
brief discussion and conclusions in Sect. [7] 


2. Observations and data reduction 

We obtained high-resolution spectra using the UVES echelle 
spectrograph (De kker et al.| |2000) mounted on the 8m UT2 
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Table 1 . Overview of observations. ’Exp. Time represents the exposure 
time. The red arm is the combination of the lower and upper part of the 
mosaic CCD chip. The parameter V represents the visual magnitude. 


Object name 

Exp. time (s) 

Exp. time (s) 

V 

(IRAC) 

Blue arm 

Red arm 

mag 

J051213.81-693537.1 

6 x 1500 

6 x 1500 

15.7 

J051848.86-700246.8 

1 x 1500 

2 x 1500 

15.5 


Kueyen Telescope of the VLT array at the Paranal Observatory 
of ESO in Chili. We selected the dichroic beam-splitter resulting 
in a wavelength coverage for the blue arm from approximately 
3280 to 4530 A, and in the red arm for the lower and upper part 
of the mosaic CCD chip from approximately 4780 to 5770 A 
and from 5800 to 6810 A, respectively. A slit width of 1 arc- 
second was used as a compromise between spectral resolution 
and slit-loss minimalisation. An overview of the observations is 
given in Table [T] For J051848.86-700246.8 (hereafter abbrevi¬ 
ated to J051848), only one blue spectrum was available because 
of an instrumental error. For J051213.81-693537.1 (hereafter ab¬ 
breviated to J051213), two additional spectra were also obtained 
but because of the significant lower quality with respect to the 
six other spectra, these two spectra are not used for this study. 
The spectra of each object were obtained in one night and in one 
large time interval. 

The data was reduced using the UVESpipeline (version 
5.3.0) in the Reflex environment of ESO [j This reduction 
includes the standard steps of extracting frames, determining 
wavelength calibration, applying this scale to flat-field divided 
data and cosmic-clipping. The standard reduction parameters of 
the UVES pipeline were used as these gave the best signal-to- 
noise (S/N) of the final spectra. 

After reduction, weighted mean spectra were calculated for 
the three wavelength ranges. Since all spectra were taken in 
small time intervals, no velocity corrections were needed for 
calculating the weighted mean. These weighted mean spectra 
were then divided into subspectra, each with a fixed wavelength 
range of 120 A for which the first and last 10 A overlapped with 
the previous and subsequent spectra. Each of these subspectra 
was normalised individually by fitting fifth order polynomials 
through interactively defined continuum points. All subspectra 
were then again merged into one large spectrum, which is used 
for the spectral analysis. In the wavelength regions where the 
subspectra overlapped, the mean flux was calculated. 


The wavelength range from 3280 to about 3900 A has a 
poor S/N, making it unsuitable for accurate spectral analysis. 
We therefore did not use this spectral region in the analysis. The 
mean S/N in the remaining wavelength range of the blue spec¬ 
trum is about 25 for both stars. The red spectra have a higher 
mean S/N between 80 and 100. The high resolution of the UVES 
spectra permits accurate radial velocity determinations using the 
positions of spectral lines. Line positions are measured by fitting 
a Gaussian through the measured line profiles, line identifica¬ 
tion is based upon the rest wavelengths from the VALD database 
|Kupka et al.| ( [1999] ). We find a radial velocity of 227.8 ±1.9 
km/s for J051848 and a radial velocity of 284.3 ±1.0 km/s for 
J051213. Both radial velocities are consistent with the average 
radial velocity of the LMC, which is approximately 270 km/s 
( [van der Marel et al.|2002| ). 


1 https ://w ww.eso.org/sci/software/reflex/ 
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Fig. 1 . Comparison of the normalised spectra of J051213 (upper), 
J051848 (middle), and J004441 (lower). The upper and lower spectra 
have been shifted in flux for clarity and all spectra are shifted to a zero 
velocity. Red and black vertical lines mark positions of 5-nuclei and non 
5-nuclei, respectively. See text for more information. 


3. Spectral analyses 


A visual inspection of the UVES spectra immediately shows the 
very strong 5-process enrichment for both objects. Figs. [T] and 
[2] show different wavelength regions of J051213 and J051848. 
Both figures also show the corresponding spectral region of 
J004441.04-732136.4 (J004441), which is a strongly 5-process 
enhanced SMC post-AGB star ( |De Smedt et al.| 2012 ). All 
three objects have similar atmospheric parameters except for the 
metallicity as described in Sect. |3.1| Both J051213 and J051848 
show strong Ba II (Fig [2]) and Y II (Fig [TJ lines, which is a 
clear indication of 5-process enrichment. Although J051213 and 
J051848 are members of the EMC, while J004441 is an SMC 
object, they have similar spectra. 

The spectral analyses includes atmospheric parameter and 
abundance determination. We use PyMOOG, our own Python 
wrapper around the local thermal equilibrium (FTE) abundance 
calculation routine MOOG (version June 2014 |Sneden|p~973| . 
For the analyses of our two EMC stars, we use the FTE Kurucz- 
Castelli atmosphere models (Castelli & Kurucz 2004[ ). As with 
the radial-velocity determination, spectral lines are identified us¬ 
ing linelists from the VAED database ( [Kupka et al.||1999| ). We 
combine the VAED line-lists with a list of lines gathered at the 
Instituut voor Sterrenkunde for the chemical analysis of A, F, 
and G stars ( [Van Winckel & Reyniers||2000| ). The line-lists im¬ 
plemented in PyMOOG cover a wavelength range from 3000 
up 11000 A. This covers the full wavelength coverage of the 
UVES spectra, and allows for the identification of spectral lines 
of about 160 ions ranging from He (Z=2) up to U (Z=92). We in¬ 
cluded neutral and firstly-ionised ions for most of the elements. 
For some 5-process elements, the second ionisation is included 
although the effective temperatures of the two sample stars in 
this study are too low for these ionisations. 

The equivalent width (EW) of spectral lines are measured in¬ 
teractively with PyMOOG. The EWs are calculated with direct 
integration. The abundances are computed with an iterative pro¬ 
cess in which the theoretical EWs of single lines are computed 
for given abundances and matched to the observational EWs. For 
our analysis, we avoided blended lines as much as possible. Syn- 
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Fig. 2. Similar to Fig. [I] except for another wavelength range. 

Table 2. Determined atmospheric parameters of J051213 and J051848. 
The errors for [Fe/H] include line-to-line scatter and model uncertainty. 
Symbols N Fe i and N Fe ii show the number of lines used for Fe I and Fe 
II, respectively. 


Object 

J051213 

J051848 

Teff (K) 

5875 ± 125 

6000 + 125 

log g (dex) 

1.00 ±0.25 

0.50 + 0.25 

ft (km/s) 

3.0 ± 0.2 

2.8 + 0.2 

[Fel/H] 

-0.56 + 0.16 

-1.06 + 0.17 

[Fell/H] 

-0.56 + 0.15 

-1.03 + 0.14 

NFei 

53 

35 

NFeii 

13 

9 


thetic spectra were used for checking whether used spectral lines 
are part of blends with other identified lines. 

3.1. Atmospheric parameters 

The atmospheric parameters are determined using the atmo¬ 
spheric parameter determination module in PyMOOG. We use 
linear interpolation to calculate atmospheric models, which are 
within the parameter steps of the Kurucz-Castelli models. For 
the determination of the effective temperature and surface grav¬ 
ity logg, we use parameter steps of 125 K and 0.5 dex, respec¬ 
tively. For the microturbulent velocity, we choose steps of 0.2 
km/s. 

To derive the atmospheric parameters, we use Fe I and Fe II 
lines. The effective temperature T e ff is determined by enforcing 
the iron abundance, derived from the individual Fe I lines, to be 
independent of lower excitation potential. We choose Fe I lines 
for this purpose since the available Fe II lines do not cover an ap¬ 
propriate range in lower excitation potential. The surface gravity 
log g is determined by enforcing ionisation equilibrium between 
the individual Fe I and Fe II abundances. The microturbulent 
velocity ft is derived by enforcing the iron abundance from indi¬ 
vidual Fe I lines to be independent of reduced equivalent width, 
which we define in this contribution as EW/T. 

The individual atmospheric parameter results together with 
the number of used individual spectral lines are shown in Table[2 
The indicated uncertainties of [Fel/H] and [Fell/H] are the tota 
errors on the iron abundances, which includes line-to-line scatter 
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Fig. 3. Spectrum synthesis of the Nd II line at 6727.695 A for J051213. 
The black spectrum is the observed spectrum, the coloured spectra rep¬ 
resent synthetic spectra with different Th abundances. The position of 
the Nd line is indicated. For more information, see text. 



Fig. 4. Spectrum synthesis of the Eu II line at 6645.064 A for J051848. 
The black spectrum is the observed spectrum, the coloured spectra rep¬ 
resent synthetic spectra with different Th abundances. The position of 
the Eu line is indicated. For more information, see text. 


and the atmospheric parameter uncertainties as described below 
in Sect. |3.2| The metallicity of J051213 is similar to the mean 
LMC metallicity of -0.4 dex, while J051848 has a significantly 
lower metallicity. 

Both stars are of spectral type F and have low surface grav¬ 
ities confirming their evolved nature. The effective temperatures 
of both stars could be determined accurately. Kamath et al. 
(2015,accepted) estimate stellar parameter for J051848, based 
on low-resolution spectra. We find that the derived stellar pa¬ 
rameters of our study and those of Kamath et al (2015, accepted) 
agree well within the estimated errors. 

3.2. Abundance determination 

We use the derived atmospheric parameters in Table [2] for the 
abundance determination using PyMOOG. We mainly use iso¬ 
lated non-blended, non-saturated lines, but this is a challenge 
because of the rich spectra with strong enhancements. However, 
with spectral synthesis, we determine the abundances of ions, 
which can only be found in blends. Lines with EWs smaller than 
5 mA are not used as they may be confused with noise in the 
spectra. 

Only a few carbon and oxygen lines are available for abun¬ 
dance determination because most of the lines are severely 
blended. For J051213, where only one spectral line was available 
for both the C and O abundances, these abundances were deter¬ 
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Fig. 5. Spectrum synthesis of the Pb I line at 4057.807 A for J051213. 
The black spectrum is the observed spectrum of J051213, the coloured 
spectra represent synthetic spectra with different Pb abundances. The 
dashed black line shows the synthetic spectrum if no Pb is present. Po¬ 
sitions of unidentified spectral lines are indicated by question marks. 
For more information, see text. 


mined using spectral synthesis. The O abundance of J051213 is 
determined using the forbidden [O I] line at 6300 A, which is 
part of a fully identified blend. This line is not sensitive to non- 
LTE effects (e.g. [Kiselman|2002| ). For J051848, we found four 
useful lines for carbon and two non-forbidden OI lines. For both 
stars, nitrogen lines larger than 5 mA are not found in the avail¬ 
able wavelength coverage of the UVES spectra. 

The spectra of both objects allowed abundance determina¬ 
tions of many elements past the iron peak. Unfortunately, all 
available Sr and Ba lines are heavily saturated hampering accu¬ 
rate abundance determinations for these two v-process elements. 
Most of the other studied v-elements abundances are determined 
from multiple, single lines. Apart from the light v-process (Is) 
peak elements Y and Zr, and heavy v-process (hs) peak elements 
La, Ce, Pr, and Nd, we find abundances for heavier elements like 
gadolinium (Gd, Z=64), dysprosium (Dy, Z=66), erbium (Er, 
Z=68), thulium (Tm, Z=69), luthetium (Lu, Z=71), and hafnium 
(Hf, Z=72). 

To illustrate our analyses, in Fig. [3] we show the compari¬ 
son between synthetic spectra with different Nd abundances for 
J051213 around the Nd II line at 6727.695 A. The red line repre¬ 
sents the best synthetic fit to the observations. The figure shows 
that our derived abundances reproduce well the observed spectral 
lines. A similar fit for the Eu II line at 6645.064 A for J051848 
is shown in Fig. [4] 

Pb is a very important tracer of the v-process nucleosynthe¬ 
sis. Unfortunately, at the photospheric conditions of our stars, the 
useful Pb lines are in the very blue part of the spectrum. Because 
of the low S/N around the strongest Pb I line at 4057.807 A, it 
is difficult to determine accurate lead (Pb, Z=82) abundances for 
both sample stars, therefore, we prefer to determine Pb abun¬ 
dance upper limits. For J051213, the S/N is poor as can be seen 
in Fig. [5] The combination of many unidentified spectral lines, 
i.e. lines not included in the line-lists, and the poor S/N make 
it difficult to determine the position of the continuum. For some 
unidentified lines, we mark the positions indicated with question 
marks. For J051848, the S/N is higher, but again many uniden¬ 
tified lines are present in the spectral region (as shown in Fig. 
[6]). For both stars, we use the spectral blended lines at 4057.5, 
4058.3, 4058.8, and 4059.9 A to estimate the position of the con¬ 
tinuum at the Pb I line by eye. With our adopted local continuum 
positions, small Pb line features are present in the observed spec- 
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Fig. 6. Spectrum synthesis of the Pb I line at 4057.807 A for J051848. 
The black spectrum is the observed spectrum, the coloured spectra rep¬ 
resent synthetic spectra with different Pb abundances. The dashed black 
line shows the synthetic spectrum if no Pb is present. Positions of 
unidentified spectral lines are indicated by question marks. For more 
information, see text. 


tra as the dashed lines represent synthetic spectra without the Pb 
I line at 4057.807 A. 

4. Abundance results 

The complete abundance analysis results of both objects are pre¬ 
sented in Tables [3] and |4] An overview of the lines used for these 
analyses can be found in two catalogues, which are available at 
CDSq] The information in each catalogue is tabulated as follows: 
Col.l lists the name of the ion to which the corresponding line 
belongs, Col. 2 shows the rest wavelength, Col. 3 contains the 
lower excitation potential, Col. 4 mentions the logarithm of the 
oscillator strength, Col. 5 lists the measured EW, and Col. 6 gives 
the deduced abundance for the corresponding line. 

Based upon the ionisation potential of the corresponding ion, 
the element over iron ratios ([X/Fe]) in Tables [3] an d0 are cal¬ 
culated using Fe I or Fe II. If the ionisation potential of an ion 
is below the ionisation potential of Fe I, the abundance of Fe 

1 is used for calculating [X/Fe]. If the ionisation potential ex¬ 
ceeds the ionisation potential of Fe I, Fe II is used for calculating 
[X/Fe]. The same principle is also used for calculating the total 
error of [X/Fe]. The errors were determined using the method 
described in [Deroo et~aL| |2005| . The uncertainties due to the 
atmospheric parameters are calculated by determining the abun¬ 
dances of a certain ion for atmospheric models with Y e ff plus and 
minus 125 K, models with log g plus and minus 0.25 dex and mi- 
croturbulent velocity plus and minus 0.2 km/s (see Table[2]). The 
uncertainty in microturbulent velocity can have important effects 
for ions, which are only determined by large lines like Y II. The 
total uncertainties are then the quadratic sum of the uncertainties 
of the mean due to line-to-line scatter (<xi 2 i), uncertainties due to 
atmospheric parameters (<xx e// , cr logg , cr^), and the Fe abundance 
uncertainty (crpe): 

= 'i(^) + + ^ +1(<rs)2+ fe) - 

The scaling to Fe I or Fe II based upon ionisation potential 
strongly decreases the [X/Fe] uncertainty of certain ions with 
respect to using only Fe I or only Fe II for all [X/Fe] calcula¬ 
tions. The [X/Fe] results in Tables [3] and [4] are plotted in Figs. [7] 
and[8j respectively. 

2 A link to the catalogues 


Table 3. Abundance results of J051213.81-693537.1. 



T e n 

log g 

= 5875 K 
= 1.00 dex 

J051213 

= 3.0 km/s 
[Fe/H] = -0.56 dex 

Sun 

Ion 

N 

log e 

n~i 2 i 

[X/Fe] 

tot 

log €q 

Cl 

1 

8.75 

0.20 

0.88 

0.26 

8.43 

OI 

1 

8.65 

0.20 

0.52 

0.26 

8.69 

Mg I 

1 

7.00 

0.20 

-0.04 

0.24 

7.60 

All 

1 

6.60 

0.20 

0.71 

0.25 

6.45 

Si I 

4 

7.32 

0.05 

0.37 

0.18 

7.51 

SI 

2 

6.51 

0.03 

-0.05 

0.13 

7.12 

Cal 

11 

6.00 

0.09 

0.22 

0.09 

6.34 

Ca II 

1 

6.00 

0.20 

0.22 

0.24 

6.34 

Sc II 

3 

2.93 

0.03 

0.34 

0.10 

3.15 

Til 

2 

4.69 

0.07 

0.30 

0.13 

4.95 

Till 

7 

4.62 

0.13 

0.23 

0.11 

4.95 

VI 

1 

3.91 

0.20 

0.54 

0.25 

3.93 

VII 

2 

3.91 

0.04 

0.54 

0.14 

3.93 

CrI 

4 

5.20 

0.07 

0.12 

0.10 

5.64 

Cr II 

7 

5.16 

0.11 

0.08 

0.10 

5.64 

Mnl 

1 

4.69 

0.20 

-0.18 

0.24 

5.43 

Fe I 

53 

6.94 

0.10 

0.00 

0.06 

7.50 

Fe II 

13 

6.94 

0.09 

0.00 

0.07 

7.50 

Nil 

4 

5.75 

0.08 

0.09 

0.12 

6.22 

Znl 

1 

4.22 

0.20 

0.22 

0.29 

4.56 

Y II 

5 

3.13 

0.11 

1.48 

0.13 

2.21 

Zr II 

2 

3.20 

0.02 

1.18 

0.10 

2.58 

Fa II 

6 

2.52 

0.13 

1.98 

0.13 

1.10 

Ce II 

6 

2.91 

0.12 

1.89 

0.13 

1.58 

Pr II 

6 

1.96 

0.06 

1.80 

0.13 

0.72 

Nd II 

14 

2.67 

0.10 

1.81 

0.13 

1.42 

Sm II 

1 

1.70 

0.20 

1.30 

0.24 

0.96 

Eu II 

2 

0.98 

0.07 

1.02 

0.13 

0.52 

Dy II 

2 

2.00 

0.20 

1.46 

0.21 

1.10 

Er II 

2 

1.90 

0.20 

1.54 

0.22 

0.92 

Pb I u 

1 

3.30 

0.20 

2.11 

0.25 

1.75 

Th II U 

1 

1.35 

0.20 

1.89 

0.27 

0.02 


Notes: The table lists for each ion the used number of lines (N) for 
the abundance determination, the determined abundance (log e) , the 
uncertainty of this abundance due to line-to-line scatter (cr m ), the el¬ 
ement over iron ratio ([X/Fe]), and total uncertainty ( cr tot ). The total 
uncertainty cr tot includes line-to-line scatter and atmospheric param¬ 
eter uncertainty (see text for details). The last column lists the solar 
abundances from |Asplund et al. 1(2009) . 

We choose a standard deviation of 0.20 dex for all ions for which the 
abundance was derived using spectrum synthesis or for which only 
one line is available. 

(u) These abundances are upper limits. 


4.1. Carbon and oxygen 

Both stars are carbon enhanced. The uncertainties of C and O 
for J051213 are large since only one line could be used for the 
abundance determinations. Both ions are also strongly sensitive 
to changes in effective temperature. We find [C/Fe] values of 
about 0.9 and 1.2 for J051213 and J051848, and [O/Fe] values 
of about 0.5 and 0.8. This results in moderate C/O ratios of only 
1.26 ± 0.40 for J051213 and 1.29 ± 0.30 for J051848, which 
are a combination of relatively low C enhancements as well as 
relatively high O enhancements. Based upon their abundances 
(Fig.0 an and their luminosities (see Sect. [5]), both stars are 
carbon-enhanced post-AGB stars but the C/O ratios are not very 
large. 
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Table 4. Abundance results for J051848.86-700246.9. 



F e ff 

log# 

= 6000 K 
= 0.50 dex 

J051848 

^ t - 2.8 km/s 
[Fe/H] = -1.03 dex 

Sun 

Ion 

N 

log e 

n~i2i 

[X/Fe] 

o-tot 

loge© 

Cl 

4 

8.61 

0.13 

1.21 

0.16 

8.43 

OI 

2 

8.50 

0.05 

0.84 

0.19 

8.69 

Mg I 

2 

6.73 

0.07 

0.19 

0.12 

7.60 

SI 

1 

6.46 

0.20 

0.37 

0.25 

7.12 

Cal 

8 

5.56 

0.08 

0.28 

0.08 

6.34 

Sc II 

6 

2.30 

0.08 

0.18 

0.11 

3.15 

Till 

4 

4.12 

0.08 

0.20 

0.11 

4.95 

CrI 

6 

4.70 

0.07 

0.12 

0.08 

5.64 

Cr II 

12 

4.78 

0.07 

0.17 

0.08 

5.64 

Fe I 

35 

6.44 

0.07 

0.00 

0.05 

7.50 

Fe II 

9 

6.47 

0.12 

0.00 

0.10 

7.50 

Nil 

8 

5.22 

0.10 

0.06 

0.10 

6.22 

Cul 

1 

3.53 

0.20 

0.40 

0.23 

4.19 

Znl 

1 

3.78 

0.20 

0.25 

0.29 

4.56 

Y II 

6 

2.79 

0.12 

1.61 

0.15 

2.21 

Zr II 

2 

2.85 

0.05 

1.30 

0.13 

2.58 

La II 

7 

2.55 

0.11 

2.48 

0.20 

1.10 

Ce II 

5 

2.37 

0.07 

1.82 

0.13 

1.58 

Pr II 

11 

1.75 

0.05 

2.06 

0.14 

0.72 

Nd II 

14 

2.69 

0.10 

2.30 

0.16 

1.42 

Sm II 

2 

1.80 

0.12 

1.87 

0.18 

0.96 

Eu II 

2 

0.81 

0.12 

1.32 

0.17 

0.52 

Gd II 

10 

1.75 

0.12 

1.71 

0.12 

1.07 

Dy II 

4 

1.91 

0.09 

1.84 

0.12 

1.10 

Er II 

1 

1.63 

0.20 

1.74 

0.26 

0.92 

Tm II 

1 

1.20 

0.20 

2.13 

0.24 

0.10 

Lu II 

3 

1.16 

0.08 

2.08 

0.15 

0.10 

Hf II 

1 

1.84 

0.20 

2.02 

0.25 

0.85 

Pb I 

1 

2.62 

0.20 

1.93 

0.25 

1.75 


Notes: Same as for Table [3] 



Fig. 7. [X/Fe] results of J051213. The errors bars represent the total 
uncertainties cr to t. The abundances of Pb (Z=82) and Th (Z=90) are 
upper limits and marked with down arrows. 


4.2. a-elements 

Concerning the available cr-elements Mg, Si, S, Ca, and Ti, 
the mean of [X/Fe] is [cr/Fe] = +0.14 and [a/ Fe] = +0.26 for 



Fig. 8. [X/Fe] results of J051848. The errors bars represent the total 
uncertainties cr to t. The abundance of Pb (Z=82) is an upper limit and 
marked with a down arrow. 


J051213 and J051848, respectively. For J051848, no Si abun¬ 
dance was available. With respect to Galactic abundances, these 
cr-element abundance are slightly deficient but they fall within 
the expected abundance ranges for the LMC, consistent with 
their respective metallicities (e.g Van der Swaelmen et al .|2013[ 
|Pompeia et al.|2008] ). 


4.3. s-process elements 

The s-process elements can be subdivided into three groups. This 
division is based upon the number of neutrons in the nuclei. The 
first group is the light s-process (Is) elements around neutron 
magic number 50 (Z ~ 38) with elements like Sr, Y, and Zr. For 
both stars, we were only able to measure Y and Zr. The second 
group is the heavy s-process (hs) elements around magic neutron 
number 82 (Z ~ 58) with elements like Ba, La, Ce, Pr, Nd, and 
Sm. The last group consists of one element, which is the double 
magic ion 208 Pb. This has a magic neutron number of 126 and a 
magic proton number of 82. Because of its double magic state, 
208 Pb is accepted as the end product of s-process nucleosynthe¬ 
sis. 

The [X/Fe] results in Figs. [7] and [8] show strong enrichment 
in both stars, confirming their post third TDU status. In particu¬ 
lar, J051848 shows strong enhancements for La (Z=57) and Nd 
(Z=60), two typical s-process elements of Ba peak. Also, the s- 
elements of the Sr peak are enriched, albeit less strongly then the 
Ba-peak elements. The abundances of elements past the hs peak 
are strongly enhanced in both stars. 

For J051213, we find a Pb abundance upper limit that is simi¬ 
lar to the overabundances of the Ba-peak elements. For J051848, 
we find that the Pb upper abundance limit is lower than the abun¬ 
dance ratios of La and Nd, and similar to the Ce (Z=58) abun¬ 
dance ratio. 

An overview of the C/O ratio, metallicity, a-element en¬ 
hancement, and the observational indices used for describing s- 
process overabundances and distributions (See Sect. [6]) are listed 
in Table [5] 
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Table 5. Overview of the C/O ratio, metallicity, ar-element enrichment, and ^-process indices for J051213 and J051848. 


Object 

C/O 

[Fe/H] 

[a/Fe] 

[ls/Fe] 

[hs/Fe] 

[s/Fe] 

[hs/ls] 

J051213 

1.26 ±0.40 

-0.56 + 0.15 

0.14 + 0.07 

1.33 ±0.08 

1.74 + 0.08 

1.61 +0.06 

0.41+0.12 

J051848 

1.29 ±0.30 

-1.03 ±0.14 

0.26 ± 0.08 

1.46 + 0.10 

2.12 + 0.08 

1.90 + 0.07 

0.66 ±0.13 


Table 6. Estimated luminosity and reddening of J051213 and J051848. 


Object 


J051213 


J051848 


E(B-V) 
L( L 0 ) 


0.69 ± 0.03 
6700 ± 200 


0.44 ± 0.02 
6250 ± 200 


5. Luminosity and initial mass determination 

5.1. Spectral energy distributions and luminosities 

The spectral energy distributions (SEDs) and known distance to 
the LMC allow us to determine the luminosities of our sample 
stars and estimate their initial masses (Sect. |5.2| ). The photomet¬ 
ric data for constructing the SEDs are retrieved from the fol¬ 
lowing catalogues: the UBVR CCD Survey of the Magellanic 
Clouds ( |Massey||2002| ), the 2MASS 6X catalogue (Cutri et al. 


2003), the Deep Near-Infrared Survey (DENIS. Fou que et al. 


(2000)), the WISE All-Sky Data catalogu e ([Cutri et al.||2012j h 


and the Spitzer SAGE survey of the LMC ( [Meixner et al.|2006] 



Fig. 9. SED of J051213. Red dots symbolise the red, original photom¬ 
etry, blue dots represent the dereddened photometry. The black line is 
the scaled Kurucz atmosphere model. 


For our sample post-AGB stars, we assume that three sources 
contribute to the reddening. The first source is reddening by 
interstellar dust in the Milky Way galaxy towards the LMC. 
[Schlegel et al.| ( |1998| find a relatively small extinction towards 
the LMC of E(B-V) = 0.075 mag. The second source is redden¬ 
ing by interstellar dust in the LMC. We assume similar extinction 
laws in the LMC and Milky Way galaxy and use the Galactic ex¬ 
tinction curves of [Cardelli et al.[fl989| ) to determine reddening 
in the LMC. The third source of reddening is caused by circum- 
stellar dust of the post-AGB object itself. We assume the cir- 
cumstellar extinction has the same wavelength dependency as 
the Galactic extinction law. 

The total dereddening is determined by applying a x 2 min- 
imalisation on the fit between the dereddened broadband fluxes 
and the appropriate Kurucz model atmospheres that are used for 
the abundance determination. For the scaled model of J051213, 
we have linearly interpolated between available models to ob¬ 
tain our preferred atmosphere model. The error on E(B-V) is 
determined by a Monte Carlo simulation of 100 arrays with a 
normal distribution of the original flux. The luminosities are cal¬ 
culated by integrating the surface of the scaled Kurucz models 
and applying a distance of 50 k pc to the LMC (e.g. Kell er &| 
Wood|2006} Reid & Parker |2010] ). The errors on the luminosity 


are determined by a Monte Carlo simulation, similar to the error 
derivation of the reddening. The estimated luminosity and red¬ 
dening results of both the sample stars are listed in Table [6] The 
constructed SEDs are shown in Figs. |9| and [T0| 

The SEDs in Figs. [9] and [TO] are double peaked, indicating 
a circumstellar dust envelope origi natin g from the strong masv- 
loss during the AGB ( [Van Winckel|2003) . We find a significant 
reddening value for both stars; in particular, J051213 has a high 
line-of-sight extinction of E(B-V) « 0.7. This extinction is domi¬ 
nated by the circumstellar extinction given the evidence for dusty 
circumstellar envelopes of these objects. The derived luminosi¬ 
ties in Table [6] are within the expected luminosity range for post- 



Fig. 10. SED of J051848. Symbols are identical to Fig. [5] 


AGB stars and they agree well with the post-TDU nature of these 
objects. 

5.2. Initial mass estimates 

Two key parameters in state-of-the-art evolution and nucle- 
osynthetic AGB models are the metallicity and initial stellar 
mass. We use the post-AGB evolutionary tracks of 
|& Wood| |1994| to obtain an estimate of the initial 
our sample stars. We compare the positions of the sample stars 
with those of the theoretical tracks in the Hertsprung-Russell 
(HR) diagram in Fig. El We have also included the posi¬ 
tions of three other v-process enriched LMC post-AGB stars: 


Vassiliadis 
masses of 
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J050632.10-714229.9, J052043.86-692341.0, and J053250.69- 


2013), and the ^-process enriched 41 


713925.8 (van Aarle et al. 

SMC post-AGB star J00444L04-732136.4 flDe Smedt et~aT 
|2012) . For all stars, we use the theoretical tracks that were calcu¬ 
lated for metallicity, which corresponded the most to the derived 
metallicities of the individual objects. The theoretical tracks 
from [Vassiliadis & Wood] |1994) start at an effective tempera¬ 
ture of 10 4 K, which is higher than the derived temperatures of 
the shown stars. Since the luminosity of post-AGB stars remains 
approximately constant when crossing the Hr-diagram, we have 
used linear extrapolation between log T e ff and log L/L© to ex¬ 
trapolate towards lower temperatures. 

We find that the ^-process enriched post-AGB stars included 
in Fig. [TT] have low masses below about 1.5 M 0 . We find that 
strong s-process enrichment is linked to low initial masses in all 
sources studied so far (|De Smedt et al.||2012| |van Aarle et al.| 
|20l3) . 

We remark that our initial mass estimates depend on the ap¬ 
plied mass-loss rate history on the AGB, and the theoretical post- 
AGB tracks of Vassili adTs & Wood] ( |1994) used for the initial 
mass estimates, are determined using initial-final mass relations 
deduced at that time. 


6. Neutron irradiation 

The ^-process distributions and s-process overabundances are 
typically represented by four observational indices: [ls/Fe], 
[hs/Fe], [s/Fe], and [hs/ls]. Unfortunately, the elemental abun¬ 
dances used for the calculations of these indices vary in the lit¬ 
erature. To have a significant statistical sample, we include the 
abundance results from |Van Winckel & Reyniers| ( |2000) and 
|Reyniers et al.| ( |2004) for Galactic objects in the ^-process in- 
dex study. Fo r the ls-index, we follow the suggestion from |Busso| 
jet al.| ( |1995) and use the mean of the relative abundances of Y 
and Zr. For the hs-index, normally the Ba abundance is used, but 
since an accurate Ba abundance study is hampered by the strong 
saturated Ba lines in ^-process enriched objects, we replace the 
Ba abundance by the Ce abundance, which can be determined 
accurately in enriched objects. The hs-index is then the mean 
of the relative abundances of La, Ce, Nd and Sm. Our four s- 
process indices are 


[ls/Fe] = 


[Y/Fe] + [Zr/Fe] 


[hs/Fe] = 


[La/Fe] + [Ce/Fe] + [Nd/Fe] + [Sm/Fe] 


[s/Fe] = 


[(ls+hs)/Fe] 
6 ; 


[hs/ls] = [hs/Fe] - [ls/Fe]. 

Since we use a non-standard element selection for the hs in¬ 
dex calculations, we list the four s-process indices of all stars in 
our s-process indices analy sis in Table [7] This includes Galac¬ 
tic object from | Van Winckel & Rey niers (2000) and Rey niers| 
et al.| (|2004), an d Magellanic Cloud objects from this paper, 


De Smedt et al.| ( |2012) , and |van Aarle et al.| ([2013 ). For ob¬ 


jects IRAS 04296+3429, IRAS 19500-1709, and J053250.69- 
713925.8, no Sm abundance could be determined. To estimate 
the Sm abundance for these objects, we scale the abundance 
of Sm to the abundance of Nd, which is the element with an 
atomic mass closest to Sm and is included in the hs index. We 
use the AGB nucleosynthesis models from the online-database 
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Fig. 11. Comparison of our sample stars and four other ^-process 
enriched Magellani c Cloud post-AGB s tars. T he evolutionary tracks 
shown are from Vassili adis & Wood] \\99A) with a metallicity of 
Z=0.001 in the upper panel and Z=0.004 in the lower panel. The mass 
corresponding to each track is indicated. For more information, see text. 


FRUITTp] (Franec Repository of Upgraded Isotopic Tables and 
Yields [Cristallo et al.1201 1) and use the models with metallicities 
closest to the respective star’s metallicity. We choose a standard 
mass of 1.5 M© based upon our estimated initial mass results for 

the Magellanic Cloud objects (see S ect.|5.2) _ 

Van Winckel & Reyniers| ( |2Q0Q) , |Reyniers et alT| ( |2004 ), and 
van Aarle et al. ( |2013) find~a strong correlation between [s/Fe] 


and [hs/ls] indices in Galactic and Magellanic Cloud objects, 
for which elements with high [s/Fe] generally show a higher 
[hs/ls] ratio. We add three objects (our two sample stars and 
J004441 from [De Smedt et al. [20l2] ) to this relation (see Fig. 
[12) . We find that these three objects confirm the correlation be¬ 
tween [s/Fe] and [hs/ls] with a correlation coefficient of 0.86. 
The [s/Fe] index represents the third dredge-up efficiency, al¬ 
though it is also influenced by the mass-loss history during the 
AGB phase and the envelope mass during TDUs. Note that the 
[hs/ls] index represents neutron irradiation, the number of neu¬ 
trons that are available for each iron seed nucleus. Generally, 
an overabundance of hs-elements with respect to ls-elements is 
expected in low-mass and low-metallicity AGB stars (1-3 M©) 


http ://fruity. oa-teramo .inaf.it/ 
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Table 7. Overview of the four ^-process indices of the stars in Figs. 12 and 13 


Object 

[Fe/H] 

C/O 

[ls/Fe] 

[hs/Fe] 

[s/Fe] 

[hs/ls] 

IRAS 04296+3429° 

-0.62 + 0.20 


1.66 + 0.12 

1.34 + 0.07 

1.45 + 0.06 

-0.32 + 0.14 

IRAS 05341+0852° 

-0.85 + 0.20 

1.4 + 0.3 

1.92 + 0.08 

2.30 + 0.07 

2.17 + 0.05 

0.38 + 0.11 

IRAS 06530-0213 fe 

-0.46 + 0.20 

2.8 + 0.3 

1.82 + 0.09 

2.11+0.07 

2.02 + 0.05 

0.29 + 0.11 

IRAS 07134+1005“ 

-1.00 + 0.20 

1.0+ 0.3 

1.55 + 0.17 

1.36 + 0.07 

1.43 + 0.07 

-0.19 + 0.18 

IRAS 08143-4406 6 

-0.39 + 0.20 

1.3 + 0.3 

1.52 + 0.15 

1.49 + 0.07 

1.50 + 0.07 

-0.03 + 0.16 

IRAS 19500-1709“ 

-0.60 + 0.20 

1.1 +0.3 

1.38 + 0.11 

0.90 + 0.07 

1.06 + 0.06 

-0.48 + 0.13 

IRAS 22223+4327“ 

-0.31 +0.20 

1.2+ 0.3 

1.55 + 0.10 

1.25 + 0.09 

1.35 + 0.07 

-0.30 + 0.14 

IRAS 23304+6147“ 

-0.79 + 0.20 

2.9 + 0.3 

1.55 + 0.10 

1.50 + 0.09 

1.52 + 0.07 

-0.05 + 0.14 

J004441.04-732136.4“ 

-1.34 + 0.17 

1.9 + 0.7 

2.06 + 0.10 

2.58 + 0.10 

2.40 + 0.10 

0.52 + 0.10 

J050632.10-714229+ 

-1.15 + 0.17 

1.5 + 0.3 

1.42 + 0.07 

1.07 + 0.14 

1.19 + 0.09 

-0.35 + 0.15 

1051213.81-693537.1“ 

-0.56 + 0.15 

1.3 + 0.4 

1.33 + 0.08 

1.74 + 0.08 

1.61 +0.06 

0.41 +0.12 

1051848.86-700246.9“ 

-1.03 + 0.14 

1.3 + 0.3 

1.46 + 0.10 

2.12 + 0.08 

1.90 + 0.07 

0.66 + 0.13 

1052043.86-692341 ,0 d 

-1.15 + 0.17 

1.6+ 0.9 

1.67 + 0.12 

1.85 + 0.10 

1.79 + 0.08 

0.19 + 0.16 

1053250.69-713925.8 rf 

-1.22 + 0.11 

2.5 + 0.7 

1.51 +0.16 

1.94 + 0.10 

1.80 + 0.08 

0.43 + 0.19 


(a) Galactic object, from 

(b) Galactic object, from 

Van Winckel & Reyniers C2000|. For C/O, we assume an uncertainty of 0.3 dex. 

Reyniers et afl(2004). For C/O, we assume an uncertainty of 0.3 dex. 

(c) SMC object, from De Smedt et al. 

(2012 

i. 

). 

(d) LMC object, from van Aarle et al. 

(2013 


(e) LMC object, from this paper. 



[s/Fe] 



[Fe/H] 


Fig. 12. Correlation between the total enrichment in ^-process elements 
and the [hs/ls] index for the stars listed in [7] Galactic objects are rep¬ 
resented by blue dots; Magellanic Cloud objects are represented by red 
diamonds. The full line shows the least-squares fit to all data points. 


Fig. 13. The absence of correlation between the metallicity [Fe/H] and 
the [hs/ls] index for all stars listed in[7] Symbols and the line are similar 
to Fig. [12] There is no clear correlation between metallicity and this 
range of metallicities. 


in which the 13 C(a,n) 16 0 is expected to be the main neutron 


source for the s-process (e.g. Straniero et al. 1995; Gallino et al. 

1998; 

Abia et al. 2002f |Cristallo et al. 20Ilf Karakas & Lat- 

tanzio 

2014). In heavier AGB stars, the ^Ne((r,n) 2 ^Mg reaction 


is dominant, which mainly produces ls-elements. Of all the Mag¬ 
ellanic Cloud objects listed in Table[7J all objects have low initial 
masses and an overabundance of hs with respect to Is, except for 
J050632.10-714229.9 (Fig. [TT). 

As described above, the [hs/ls] index is expected to increase 
with decreasing metallicity because of a larger number of neu¬ 
trons for each iron seed, assuming that the diffusion of protons 
in the He-rich intershell is totally independent of metallicity. In 
Fig.[l3] which also includes our sample stars and J004441, there 
is no clear correlation between [Fe/H] and [hs/ls]. 

Low-mass, low-metallicity stars are expected to have high 
overabundances of Pb with respect to the hs-elements. The AGB 
models predict the 13 C neutron source to be of primary origin, 
which means that it is independent of the initial metallicity. For 


similar models with different metallicity, the neutron production 
is predicted to be similar, but the number of available Fe seeds 
is different. Per Fe seed nucleus, there are more neutrons avail¬ 
able at low metallicity. This results in the creation of heavier s- 
elements. Since Pb is the expected end-product of the v-process 
nucleosynthesis, high Pb abundances are expected in metal-poor 
stars. 

In order to compare our Pb abundance results to mean the¬ 
oretical predictions, we obtained predicted values from differ¬ 
ent AGB evolution and nucleosynthesis codes. For this compar¬ 
ison, we use [Pb/hs] ([Pb/Fe] - [hs/Fe]) and [Pb/ls] ([Pb/Fe] - 
[ls/Fe]), which represent the overabundance of Pb with respect 
to the other v-process elements. Fig. [14] shows the comparison 
between the observed [Pb/hs] (upper panel) and [Pb/ls] (lower) 
upper limits, together with 2.0 M© model predictions of the 
STAREVOL code ( |Siess|2007] and references therein) for metal¬ 
licities of [Fe/H] = -1.0, -0.5 and -0.2 dex (red fu ll line), Mount- 
Stromlo Evolutionary code (MSE) predictions ( [Fishlock et al.[ 
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Fig. 14. Overview of the observed [Pb/hs] (upper panel) and [Pb/ls] 
(lower panel) upper limits of ^-process enriched post-AGB stars in the 
Magellanic Clouds from this study and |De Smedt et al.] ( |2014| ). The ob¬ 
served abundance upper limits are plotted together with the [Pb/hs] and 
[Pb/ls] predictions of the 2.0 M 0 STAREVOL models (red full line),the 
2.0 M 0 and 2.1 M 0 Mount Stromlo models (black crosses and yellow 
diamond, respectively) and the 1.5 and 2.0 M 0 FRUITY models (green 
dashed and blue dot-dashed lines, respectively). The black horizontal 
lines represent the [Fe/H] uncertainty of the stars shown. The position 
of each star is indicated with the first part of its 2MASS name. For more 
information, see text. 


12014| [Karakas|2010| and references therein) for [Fe/H] « -1.3 
and -0.15 dex for 2.0 M© models (black crosses) and [Fe/H] « 
-0.5 dex for a 2.1 M© model (yellow diamond), and the publicly- 
available 2.0 M© model predictions of FRUITY ( |Cristallo et al.| 
2011) within a metallicity range from [Fe/H] = -1.5 up to 0.0 


dex (blue point-dashed line) for seven different metallicities. Be¬ 
cause of the limited number of available metallicities within the 
shown metallicity ranges, the STAREVOL and FRUITY curves 
in Fig. [14] are not smooth. We remark that, although all of the 
stars in Fig. [l4]have estimated i nitial masses between 1.0 and 1.5 
M© (see|De Smedt et al.|2012[ van Aarle et al.|2013| and Sect. 

it is justified to make a comparison with 2.0 M© model 
predictions since the differences between the 1.5 and 2.0 Mo 
FRUITY model predictions is only marginal. 

Fig. [14] shows that the predictions fit well the observed Pb 
abundance upper limit of J051213 but strongly overestimate the 
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Pb overabundances of the lower metallicity stars. There is an in¬ 
creasing discrepancy between predicted and observed Pb abun¬ 
dances towards lower metallicities. A larger sample of post-AGB 
stars with derived Pb abundances is required to confirm the ob¬ 
served discrepancy towards lower metallicities. More research is 
clearly needed to explain the consistent low Pb abundances in v- 
process enriched low-metallicity post-AGB stars. This is beyond 
the scope of this study. 


7. Conclusions 

We aim to understand the chemical diversity demonstrated by 
post-AGB stars, and we focus on v-process rich objects. With de¬ 
tailed studies of a well-sampled grid of post-AGB objects cov¬ 
ering a wide range in luminosity (hence core-mass and initial 
mass) and metallicity, we aim to build up a systematic set of 
constraints on AGB models. These will help us to understand 
the physical processes related to AGB nucleosynthesis and mix¬ 
ing. To obtain these constraints, we study newly discovered post- 
AGB stars in the Magellanic clouds. Because of their known dis¬ 
tances, accurate luminosity and, hence, initial mass estimates can 
be made, which is contrary to the situation for Galactic peers. 

Here, we report on a detailed chemical study of two v- 
process enriched optically visible post-AGB stars in the LMC. 
The objects were carefully selected from our low-resolution sur¬ 
veys ( [van Aarle et al.|20fT]|Kamath et al.|20l5] accepted). Our 
high-resolution spectral data reveal very rich spectra, literally 
swamped with atomic lines of v-process elements. We quanti¬ 
fied that both objects belong to the most ^-process rich objects 
known to date despite their similarly low C/O ratios. The metal¬ 
licity difference between both stars is about 0.5 dex. 

With the addition of our two sample stars, all v-process en¬ 
hanced post-AGB stars in the LMC and SMC studied thus far 
cluster in the same region of the HR-diagram. All are associated 
with the latest evolutionary phase of low-mass stars with, on av¬ 
erage, low metallicity. 

We also confirm the correlation between [hs/ls] and [s/Fe] 
for v-process enriched post-AGB stars and find that both Galac¬ 
tic and Magellanic Cloud objects follow the same relation. The 
neutron exposure as traced by [hs/ls] is correlated with the over¬ 
all v-process overabundances, as traced by [s/Fe]. [s/Fe] is also 
determined by the amount of dredged-up intershell material and 
the envelope dilution. In addition, we strengthen and confirm 
that there is no clear correlation between neutron irradiation and 
metallicity. Furthermore, we find an increasing discrepancy be¬ 
tween observed and predicted Pb overabundances towards lower 
metallicities: the higher metallicity J051213 with [Fe/H] « -0.56 
dex fits well the model predictions, while J051848 and other 
Magellanic Cloud objects from previous studies reveal much 
lower [Pb/hs] and [Pb/ls] results than predicted. 

Our current and future research involves expanding our sur¬ 
vey of chemical studies in the Magellanic clouds to construct 
a grid that cover a larger and well-sampled range in luminosity 
(hence initial mass) and metallicity. We intend to use this grid of 
constraints to fine-tune AGB models and therefore the underly¬ 
ing physics of nucleosynthesis and mixing. 
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